STABILITY OF THE PERIODIC TODA LATTICE IN THE 
SOLITON REGION 



HELGE KRUGER AND GERALD TESCHL 



Abstract. We apply the method of nonUnear steepest descent to compute 
the long-time asymptotics of the periodic (and slightly more generally of the 
quasi-periodic finite-gap) Toda lattice for decaying initial data in the soliton 
region. In addition, we show how to reduce the problem in the remaining 
region to the known case without solitons. 



1. Introduction 

Consider the doubly infinite Toda lattice in Flaschka's variables (see e.g. [21] 
[221, or [21) 



b{n,t)^2{a{n,tf -a{n-l,tf), 
d{n, t) — ain, t)(b{n + l,t) — b{n, t)), 

{n, e Z X K, where the dot denotes differentiation with respect to time. We will 
consider a quasi-periodic finite-gap background solution (og, bq), to be described in 
the next section, plus a short range perturbation (a, b) satisfying 

(1.2) ^(1 -t- \n\f+\\a{n,t) - ag{n,t)\ + \b{n,t) - bg{n,t)\) < oo. 

n 

for some I E N. It suffices to check this condition for one t e R (see |6]). The 
perturbed solution can be computed via the inverse scattering transform. The case 
where {aq,bq) is constant is classical (see again [21] or [21]) and the more general 
case we want to apply here was solved only recently in [B] (see also [TS]). The 
long-time asymptotics in the case where {aq,bq) is constant were first computed 
by Novokshenov and Habibullin ^19j and later made rigorous by Kamvissis |12j 
under the additional assumption that no solitons are present. The case of solitons 
was recently investigated by us in |16j . For a self-contained introduction, includ- 
ing further references and a more detailed history of this problem, see our review 
[17] ■ The long-time asymptotic in the present situation were first established by 
Kamvissis and Teschl in [T4^ (see also jT3] for a short overview) by a generalization 
of the so-called nonlinear stationary phase/steepest descent method for oscillatory 
Riemann-Hilbert problem deformations to Riemann surfaces. While [14] contains 
only leading order asymptotics, higher asymptotics were given in [15] . However, 
both [14| and [15] assume that no solitons are present and hence the purpose of the 
present paper is to show how the results can be extended to cover solitons as well. 
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To fix our background solution, choose a Riemann surface M as in ()2.ip . a 
Dirichlet divisor T>fi, and introduce 



(1.3) z{n,t) = Ae^{^+) - aE,{'Dp.) - nA^_ (00+) + tlL, - E^^ E 



where A^^ fefio) Abel's map (for divisors) and S^^, f/g are some constants 
defined in Section [2] Then our background solution is given in terms of Riemann 
theta functions by 



2 _ .20izin + l,t))9iz{n~l,t)) 



aq{n, t) = a 



6{z{n,t)Y 



(1.4) b^(^n,t) = b+--\o^' ^-^ >' 



2dt ^\e{z{n-l,t))J' 



where a, b are again some constants. We remark that this class contains all periodic 
solutions as a special case. 

In order to state our main result, we begin by recalling that the sequences a(n, t), 
b{n, t), n € for fixed < S M, are uniquely determined by their scattering data, that 
is, by the right reflection coefficient R+{X,t), A € a{Hg), and the eigenvalues pj G 
M\a{Hq), j = 1, . . . , A'^, together with the corresponding right norming constants 
7+j(t) > 0, j — 1,...,N. Here cr{Hq) denotes the finite-band spectrum of the 
underlying Lax operator Hq. 

The relation between the energy A of the underlying Lax operator Hq and the 
propagation speed at which the corresponding parts of the Toda lattice travel is 
given by 



(1.5) v{X) ^ 



where 



(1.6) i;(A) = lim 



and can be regarded as a nonlinear analog of the classical dispersion relation. Here 
00 _ is an Abelian differential of the third kind defined in (|2.14p and f2o is an 
Abelian differential of the second kind defined in (|2.15p . We will show in Section O 
that w is a homeomorphism of R and we will denote its inverse by C,{n/t). 
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Define the following limiting lattice 
(1.7) 

' ai{j,t) y 0{z{n,t)) e{z(n-l,t) + 5{n,t)) 
,aq{j,t)) ^eiz{n^l,t)) eiz{n,t) + 6in,t)) 



n 



J =71 



Y]^ exp ( -2 / uJoo+ oo_ 

X exp ( -^T / log(l - \R\^)uoo+ oo_ I , 
\27ri 7c'(„/t) / 



E / ^0- 

1 d 



2 



)0 



2ds '^V dUin^s)) 

6,in/t)^2 E AE,iPk) + ^l log(l-|i?|^ 

Pfc<C(«/t) ^ ' ' 

where R = R+{X, t) is the associated reflection coefficient, d is a canonical basis of 
holomorphic diS'erentials, and C{n/t) — TT~^{a{Hq) n (— oo, oriented such 

that the upper sheet is to the left. 

Theorem 1.1. Assume (|1.2p and abbreviate by — v{pif.) the velocity of the k'th 
soliton defined above. Then the asymptotics in the soliton region, {{n,t)\C,{n/t) € 
'E\a{Hq)} , are as follows. 

Let e > sufficiently small such that the intervals [c^ ~ e,Ck + s], 1 < k < N, 
are disjoint and lie inside v{M.\a{Hq)). 

If \j — Ck\ < e for some k, the solution is asymptotically given by a one-soliton 
solution on top of the limiting lattice: 

n 44 - ( l ^^^^^i^Mp^^-~]f + or')] , 

t^ Q\ X^hC+\ hf^\ t ^■, aiin,t)i:i{pk,n,t)i:i{pk,n + l,t) 

(1.8) 2^ b{j,t) -bi{j,t) ^ -jk{n,t) +0{t ), 

for any I > 1, where 

oc 

(1-9) Q,T,(p,n,t) = 1 +7 E i'l.+ ipJ^tf 

j=n+l 

and 

(1.10) 7k{n,t)^^,^ffil^ 

T(Pk,n,t) 

Here ipi(j),n,t) is the Baker-Akhiezer function (cf. Section\^ corresponding to the 
limiting lattice defined above. 
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If \ j — Ck\ > £, for all k, the solution is asymptotically close to the limiting 
lattice: 

OQ 

(1-11) E h{],t)-bi{j,t)=0{t-% 

for any I >1. 

In particular, we see that the solution splits into a sum of independent solitons 
where the presence of the other solitons and the radiation part corresponding to 
the continuous spectrum manifests itself in phase shifts given by (jl.lOp . While in 
the constant background case this result is classical, we are not aware of proof even 
in the case of a pure soliton solution. Moreover, observe that in the periodic case 
considered here one can have a stationary soliton (see the discussion in Section [SJ . 

The proof will be given at the end of Section [T] Furthermore, in the remaining 
regions the analysis in Section [7] also shows that the Riemann-Hilbert problem re- 
duces to one without solitons. In fact, away from the soliton region, the asymptotics 
can be computed as in [15j . The only difference being, that in the final answer, the 
limiting lattice has to be replaced by the one defined here and the Blaschke fac- 
tors corresponding to the eigenvalues have to be added to the partial transmission 
coefficient. 

Finally, we note that the same proof works even if there are different spatial 
asymptotics as n ^ ±cx) as long as they lie in the same isospectral class (cf. 
Remark 13.31 below) . 

2. Algebro-geometric quasi-periodic finite-gap solutions 

As in [14], we state some facts on our background solution (a^, hq) which we want 
to choose from the class of algebro-geometric quasi-periodic finite-gap solutions, 
that is the class of stationary solutions of the Toda hierarchy, 2 . In particular, 
this class contains all periodic solutions. Wc will use the same notation as in [3T], 
where we also refer to for proofs. As a reference for Riemann surfaces in this context 
we recommend [9,. 

To set the stage let M be the Riemann surface associated with the following 
function 

23+1 

(2.1) Rlil^i^), R2g+2{Z) = n - ^j)' Ea<Ei<---< E2g+U 

j=o 

g G N. M is a compact, hyperelliptic Riemann surface of genus g. We will choose 
i?2g+2(^) the fixed branch 

2g+l 

(2.2) <+2(^) = - n v'^^' 

where is the standard root with branch cut along (— cx3,0). 

1 /2 

A point on M is denoted hy p = {z, ±i?2g+2(^)) — z & C, or p = (oo, ±) — 

oo±, and the projection onto CU {oo} by 7r(p) = z. The points {{Ej, 0),0 < j < 
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25 + 1} C M are called branch points and the sets 

(2.3) n± = {{z, ±Rygl^{z)) I ^ e C \ y [E2i,E2i+^]} c M 

are called upper, lower sheet, respectively. 

Let {aj, &j}j=i be loops on the surface M representing the canonical generators 
of the fundamental group 7ri(M). Wc require aj to surround the points i?2j-i, E2j 
(thereby changing sheets twice) and bj to surround -Eq) E^j-i counterclockwise on 
the upper sheet, with pairwise intersection indices given by 

(2.4) Qioaj = biobj = 0, aiobj=5ij, l<i,j<g- 

The corresponding canonical basis {Cj}j=i for the space of holomorphic differentials 
can be constructed by 

(2-5) C = E^0')^I7^' 

where the constants c(.) are given by 



Cj (fc) — Cj^ , Cjk — I 

J a 



pl/2 / jjl/2 I ■ 

The differentials fulfill 

(2.6) Ck = Sj^k, Ck= rj,k, Tj^k = Tk,j, ^ <j,k<g. 

J aj J bj 

Now pick g numbers (the Dirichlet eigenvalues) 

(2.7) (A,)f=i = (M„^,)?=i 

whose projections lie in the spectral gaps, that is, G [E2j-i,E2j]. Associated 
with these numbers is the divisor which is one at the points fij and zero else. 
Using this divisor we introduce 

z{p, n, t) = (p) - a^o (^a) - nA^_ (00+) + tUj^ - E^^ G 

(2.8) z{n,t) = z{oo+,n,t), 

where is the vector of Riemann constants 

(^•y) -Eo,j ^ , 

UjQ are the 6-periods of the Abelian differential fio defined below, and A^^^ {&Eo ) 
Abel's map (for divisors). The hat indicates that we regard it as a (single- valued) 
map from M (the fundamental polygon associated with M by cutting along the a 
and b cycles) to C^. We recall that the function 0{z{p,n,t)) has precisely g zeros 
fij{n,t) (with fij (0,0) = fij), where 

(2.10) 0{z) = exp27ri ({m,z) + ^-'^^A , 2 € C^, 

is the Riemann theta function associated with M. 
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Then our background solution is given by 
aq{n, tY = a 



e{z{n,t)Y 

The constants a, h depend only on the Riemann surface (see [21] Section 9.2]). 
Introduce the time dependent Baker- Akhiezer function 



(2.12) V,(P, t) = C{n, 0, t) fY' exp (n / u:^^oo^+t f^o), 
where C{n, 0, i) is real-valued, 

(2.13) CMtr- %(0-0))^U(-l'0)) 



6i(z(n,t))6l(z(n- l,t))' 
and the sign has to be chosen in accordance with aq{n, t). Here 

(2.14) u;oo+ oo_ = -di, 

is the Abelian differential of the third kind with poles at cxd_|_ and cxd_ and 

(2.15) ^o^nk^^rf., E^.4£'^- 

-"^23-1-2 j=0 j=0 

is the Abelian differential of the second kind with second order poles at oo+ and 
oo_ (see Sects. 13.1, 13.2]). All Abelian differentials are normalized to have 
vanishing aj periods. 

We will also need the Blaschke factor 

where E{p) is Eq if p < Eq, either i?2j-i or i?2j if p G (i?2j-i, -B2j), 1 < j < 5, and 
E2g+i if p > E-ig+i- It is a multivalued function with a simple zero at p and simple 
pole at p* satisfying \B{p,p)\ = 1, p g 911+. It is real-valued for 7r(p) g (— oo,i?o) 
and satisfies 

(2.17) B{Eo,p)^l and B{p\ p) = B{p, p*) = B{p, p)-^ 
(see e.g., [SS]). 

The Baker-Akhiezer function is a meromorphic function on M \ {oo±} with an 
essential singularity at oo±. The two branches are denoted by 

(2.18) tpq^±{z,n,t) ^ ipg{p,n,t), p^{z,±) 
and it satisfies 

Hq{t)tpq{p,n,t) = 7r(p)V'q(p,n,i), 

(2.19) ^^%{p,n,t) ^ Pq^2{t)4'qip,n,t), 
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where 

Hq{t)'ilj{n) = aq{n, t)^j{n + 1) + aq{n — 1, t)il;{n - 1) + bq{n)xj}{n), 

(2.20) Pq,2{t)i'{n) = a,(n, t)^{n + 1) - a,(n - 1, i)V'(n - 1), 
are the operators from the Lax pair, 

(2.21) j^HS) = Hq{t)Pq^,{t) - Pq,2{t)Hq{t), 

for the Toda lattice. 

It is wch known that the spectrum of Hq(t) is time independent and consists of 
g + I bands 

9 

(2.22) aiHq)=\J[Ey,Ey+i]. 

3=0 

For further information and proofs we refer to PT, Chap. 9 and Sect. 13.2]. 
3. The Inverse scattering transform and the Riemann-Hilbert 

PROBLEM 

In this section our notation and results are taken from [S] and [5] . Let '09, ± (z, n, t) 
be the branches of the Baker- Akhiezer function defined in the previous section. Let 
ip±{z,n,t) be the Jost functions for the perturbed problem defined by 

(3.1) lim w{z)^"'{ijj±{z,n,t) -ipq±{z,n,t)) ^0, 

n — '±00 ' 

where w{z) is the quasimomentum map 

(3.2) w{z) = exp{ Woo+oo_), p={z,+). 

J Eq 

The asymptotics of the two projections of the Jost function are 
ilj±{z,n,t) = , f X 

n 

(3.3) X (1+ (i?±(n,t)±^6,(j- ;,i))- + 0(-2) 
as 2: — > 00, where 

n— 1 / . ,\ n — 1 

A-M^Yl ^f-j- B.{n,t)= {bqU,t)-b{j,t)). 



(3.4) 



One has the scattering relations 



(3.5) T{z)^^{z,n,t) = 'il}±{z,n,t) + R±{z)^±{z,n,t), z & a{Hq), 

where T{z), R±{z) are the transmission respectively reflection coefficients. Here 
'4'±{z, n, t) is defined such that ip±{z, n, t) = lim^jo ■*/'±(^ + is? ^ ^ <^{Hq)- If 
we take the limit from the other side we have ^'±(•2, n, t) = lim£|o '0±(2 — is, t). 

The transmission and reflection coefficients have the following well-known prop- 
erties [5]: 
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Lemma 3.1. The transmission coefficient T{z) has a meromorphic extension to 
C\cr{Hq) with simple poles at the eigenvalues pj. The residues ofT{z) are given by 



(3.6) Resp^ T{z) 



-^2g+2(Pj) 7±J 
Pj ^ j -i-ig I \ ±1 ' 



where 

{3-7) 7i,',=ElV'±(P.-">*)P 

and ip- {pj ,n,t) — Cj ip+ {pj ,n,t). 
Moreover, 



(3.8) T{z)R+{z) + T{z)R-{z) = 0, iTiz)^ + \R±{z)f ^ 1. 

In particular one reflection coeflicient, say R{z) — R^{z), and one set of norming 
constants, say 7^ = 7+j , suffices. 

We win define a Riemann-Hilbert problem on the Riemann surface M as follows: 

We are interested in the jump condition of m{p,n,t) on S, the boundary of II-i- 
(oriented counterclockwise when viewed from top sheet n_|_). It consists of two 
copies I]± of <j{Hq) which correspond to non-tangential limits from p = (z, +) with 
±Im(z) > 0, respectively to non-tangential limits fromp = (z, — ) with =Flm(z) > 0. 

To formulate our jump condition we use the following convention: When repre- 
senting functions on S, the lower subscript denotes the non-tangential limit from 
n+ or n_, respectively, 

(3.10) rn±{pa) = lim m(p), po G S. 

Using the notation above implicitly assumes that these limits exist in the sense that 
m{jp) extends to a continuous function on the boundary away from the band edges. 

Moreover, we will also use symmetries with respect to the the sheet exchange 
map 



(3.11) P* = 

and complex conjugation 



(z,T) forp=(z,±), 
ooip for p = oo±, 



(3.12) p 



'(z,±) forp=(z,±)^S, 
(z, t) for p = (z, ±) e S, 
oo± for p = oo±. 

In particular, we have p — p* for p G E. 

Note that we have m±{p) — m^{p*) for rh{p) — m{p*) (since * reverses the 
orientation of E) and rri±{p) = rn±{p*) for ■m{p) — m{p). 

With this notation, using (13. 5p and (13. 8p . we obtain 

m+{p, n, t) — m_ (p, n, t)J{p, n, t) 



ni3) J(vnt)-( -i?(p)e(p,n,t)e-*^W\ 

(3.13) -^^P'"'^^- l^i?(p)e(p,„,i)e**(P) 1 ) 
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where 



and 



e{z{p,n,t)) e{z{p\Q,{))) 



(3.14) (/.(p,-) = 2 / no + 2- 

foT p€ S. Note 



t -'eo i JEa 



Vq[P*,n, t) 

Here we have extend our definition of T to S such that it is equal to T{z) on 
and equal to T(z) on Similarly for R{z). In particular, the condition on S+ 
is just the complex conjugate of the one on S_ since we have R{p*) — R{p) and 
m±{p*,n, t) — m±{p, n, t) for p G S. 

Furthermore, 
(3.15) 

m{p,n,t)=(A+(n,t){l-2B+(n-l,t)\) -^-i-^(l + 2i?+(n, <)i)) + 0(1), 

for p — (z, +) — s- cxD^, with A±[n^t) and B±{n,t) are defined in (|3.4p . The formula 
near cxd_ follows by flipping the columns. Here we have used 

(3.16) nz) = AAn,t)A^{n^t){l - + h{n^t) + B.{n,t) ^ 1 

z z'^ 

Using the properties of ip±{z,n,t) and ipq^±{z,n^t) one checks that its divisor sat- 
isfies 

(3.17) (mi) > -X'/i(„,t). - Pp, (m2) > -X'^Cn,*) - Vp*, 
where 

(3.18) I?p = ^I?p^, P,. 

3 3 

Here (/) denotes the divisor of /. 

Theorem 3.2 (Vector Riemann-Hilbert problem). Let S+{H{Q)) = {i?(A), A S 
(7{Hq); {pj,"fj), 1 < i < -/V} t/ie ri^/if scattering data of the operator H{0). Then 
m{z) — m{z,n,t) defined in (j3.9[) is meromorphic away from S and satisfies: 

(i) T/ie jump condition 

(3.19) 

/ ^ 7^ ^ -i?(p)e(p,n,t)e-**(f)\ 

/or p e S, 

(ii) the divisor 

(3.20) (mi) > -I?^(„,t). - Pp, (ma) > -I?p(„,t) - I?p. 
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and pole conditions 

(3.21) 

1 /2 

+ TTs— 7 rT~* 7T'^2b > -X>/i(„,t)-, near p^, 

(iii) the symmetry condition 

(3.22) m(p*) = m(p) (^J J 

(iv) and i/ie normalization 

(3.23) mi(oo+) • m2(oo-|_) = 1 mi(cx3+) > 0. 

Proof. We already derived the jump condition (|3.19p . The pole conditions follow 
since T{p) is meromorphic in M\E with simple poles at pj and residues given by 
p.6p . The symmetry condition holds by construction and the normalization p.23p 
is immediate from p.lSp . □ 

Remark 3.3. We note that the same proof works even if there are different spatial 
asymptotics as n ^ ±00 as long as they lie in the same isospectral class. In 
fact, following [7], we now have two different background operators and we will 
denote the corresponding Baker-Akhiezer functions by ^p^{p,n,t). Using 

in place of (|3.9p one easily checks that Theorem \3.2\ still holds in this case. The 
only difference now is that the (right) scattering data pj and R+{z) will not satisfy 
the algebraic constraints given in 23J. 

For our further analysis it will be convenient to rewrite the pole condition as a 
jump condition. Choose e so small that the discs \'k{p) — pj\ < e are inside the upper 
sheet and do not intersect. Then redefine m in a neighborhood of pj respectively 
p* according to 



(3.25) m{p) 



1 0' 



k(p)-P3l<e 



ij ■^gip,n,t) ^ \ ' pen 



m(p) ^ 

nfc = i(Pj-Mfe) T(p)-Pj i'gip',n,t) 

I 1 ^2g + 2(P3) 7j ■4'q(P,n,t) . 

1 



^m{p), else. 
Then a straightforward calculation shows 

Lemma 3.4. Suppose m{p) is redefined as in p.25p . Then m{z) is meromor- 
phic away from E and satisfies p.l9p . (j3.22p . (j3.23p . the divisor condition change 
according to 

(3.26) {mi) > ~Vf,i^n,t)', (W2) > -X'^(„,t) 
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and the pole conditions are replaced by the jump conditions 



(3.27) 



/ 1 0\ 

VOfc^ife-Mfe) ^{P)-Pj 4>g(p',n,t) ) 

ni+{p) = m-(p) Yll^iiPj-p-k) ■^(p)-pj ■i'qip' ■n,t) I p e ^^(p*), 

Vo 1 y ' 



where 

(3.28) E,(p) = {q € n± : |7r(g) - z| = e}, p = (z, ±), 

is a small circle around p on the same sheet as p. It is oriented counterclockwise 
on the upper sheet and clockwise on the lower sheet. 

Next we turn to uniqueness of the solution of this vector Riemann-Hilbert prob- 
lem. This will also explain the reason for our symmetry condition. We begin by 
observing that if there are 5 + 1 points pj £ M, such that m(pj) — {0 0) , then by 
Riemann-Roch ([ST, Thm. A.2]), 

g+l g+l g+l 

(3.29) r{- ^ V,^ ) = deg(^ P.J + 1 - 5 + z(^ P^J > 2, 

there are at least two linearly independent functions h{p) with (h) > — X^j-^Pj- 
In particular, there is a non-constant function h{p) and n(p) = h{p)m(p) satisfies 
the same jump and pole conditions as m{p). However, it will in general violate 
the symmetry condition! In fact, the symmetry condition (|3.22p requires h{p) = 
h{p*). Hence h{p) = h{TT{p)) is the lift of a meromorphic function h{z) on C. In 
particular, ii pj is a branch point, it will have an even order pole. So if we choose 
the points Pj to be different branch points, this yields a contradiction since at least 
one of the points pj must be a non-removable pole. Thus, without the symmetry 
condition, the solution of our vector Riemann-Hilbert problem will not be unique 
in such a situation. Moreover, such a situation can be indeed created starting with 
(e.g.) a one soliton solution and using the Dirichlct commutation method ([ID]) [HI 
Sect. 11.8]) to place all 5 -I- 1 Dirichlet eigenvalues at the band edges (which equal 
the branch points of our Riemann surface) . 

Lemma 3.5 (One soliton solution). Suppose there is only one eigenvalue and a 
vanishing reflection coefficient, that is, S^{H(t)) = {R{p) = 0, p G S; {p,"f)}. Let 

00 

(3.30) Cg^y{p,n,t) = 1+7 X! ■^'/(^'i'*)^ = ^ + 'lWq,{n,t){lp'q{P, ■,t),1pq{p, ■,t)) 



j=n+l 



and 
(3.31) 

V'g,^(p,n,t) 



y{p,n,t)'lpq{p, n,t) + j^lljq{p,n, i)W^g,(n,t)(V'g(P, ■,'t),1pq{p, .,t)) 

\Jcqa{P-, n-\, t)cq,^{p, n, t) 



p = (z,±). Here Wq^(^n,t){f, o) = aq{n,t){f{n)g{n + 1) - f{n + l)g{n)) is the usual 
Wronski determinant and the prime denotes a derivate with respect to p. 
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Then the unique solution of the Riemann-Hilbert problem (13.19^ - (|3.23l) is given 

by 

mo{p)^{fip*,n,t) fip,n,t)), /(p, t) = 

In particular, 
(3.32) 

Cq.-^{p,n- l,t) ag{n,t)il;g{p,n,t)il;q{p,n+l,t) 

A+(n,t) = A -r— , B+[n,t) ^ -J . 

y Cq.f[p,n,t) 2cq,j{p,n,t) 

Proof. This follows directly from p.9p after inserting the formulas for the Jost 
functions obtained from (e.g.) the double commutation method found in [TD] or 
[m Sect. 14.5] (cf. also [8]). 

Alternatively, it can also be easily checked directly, that p.32p solves p.l9p - 
()3.23|) . In fact, except for the pole conditions everything is straightforward. The 
pole conditions follow from 

lim(z - p)fip*) = ^^^' "' ^^ WMp, t), V',(P*, i)) 

P^P \/Cq,'y{p, 71-1, t)Cg^y(p, n, t) 

1/2 



where 



and 



^ l{P-,n,t) ^2g+2(p) 

7(p, n, t) = 7T#^ = 7e(p, n, t)e'^(rf , 
il^qip* ,n,t) 



lim /(p) 



Cg.7(P, t) + iWnii'qip, .,t),Tp'g{p, .,t)) 



P^P ^Cq^~^{p,n~ l,t)Cq^^f{p,n,t) 

1 



y/cq,^{p, n - 1, t)cq^^{p, n, t) 

The formulas for A_|_(n, t) and B^{n, t) follow after expanding around p — 00+ and 
comparing with (|3TT5| . 

To see uniqueness, let rho{p) be a second solution which must be of the form 
™o(p) = if{p*) fip)) by the symmetry condition. Since the divisor I?/i(„^t)* + 'Dp 
is nonspecial ([21] Lem. A. 20]), that is, its index of speciality vanishes, i{T>jx(j^ ty + 
Vp) = 0. Thus the Riemann-Roch theorem ([211 Thm. A.2]), 

^(-^'M",*)* ~ ^p) = deg(X»^(„,t). + Pp) + 1 - .g + i{Vfj^n,tr + ^p) = 2, 
implies that there are two linearly independent functions satisfying (/) > — I'^(„jt)* — 
Vp. One is f{p) and the other one is the constant function. This implies f{p) = 
'^f{p) + P fo'' some a,/3 € C. But the pole condition implies /3 = and the 
normalization condition implies a — I. □ 

Note 



f , ,x \J Cq.'yiP, n-1, t)Cq.^{p, 71 + 1, t) 

aq^^[n,t) = aq[n,t) 



(3.33) 6,,-^(7i,t) = 6q(?i,i) +79' 



Cq,'rip,n,t) 

aq{n, t)lpq{p, 71, t)-0g(p, 71 + 1, ^) 
Cq,'y{p,n,t) 



STABILITY OF THE PERIODIC TODA LATTICE 



13 



where d*u{n) — u{n — 1) — u{n), and that aq^y(n,t), hq^-y{n^t) is centered at 

(3.34) 2a{p) {n - v{p)t) + ln{-f) = 0, 
where 

(3.35) a{p) = Re / Woo+oo^, v{p) = r^^e / ilo- 

Jeo "(p) Jeo 

Note, however, that i?<j,7 looks asymptotically like Hq as n — > +oo but not as 
-oo (cf. [8]). 

Since / has g+1 poles, there are also 5 + 1 zeros which are given by the Dirichlet 
eigenvalues of Hq^{t). Hence there is precisely one in the closure of each interior 
gap of the spectrum a{Hq^^) — cr{Hq) U {p}. Moreover, observe f{pi) = f{Pi) — 
if and only if Wn{ipq.^+{p),ipq.,+{zi)) = W„(-0g,+ (p), "0,,- (zi)) = 0. Hence we also 
must have Wn{^pq,+ {zi),ipq,-izi)) = 0, that is, zi G {Ej}. Furthermore, even in the 
general case m(pi) — (O O) can only occur at zi G {Ej} as the following lemma 
shows. 

Lemma 3.6. For m{p) defined in (|3.9p set 

If rh{pi) = (0 0) for m defined as in (j3.9p . then z\ G \Ej\- Moreover, the zero 
of at least one component is simple, in the sense that mk{p)^^ = 0{{z — zi)^^^^), 
in this case. 

Note that the theta functions in (|3.36p are used to cancel the poles at fi{n,t) 
respectively jl{n,t)* such that their contribution is separated from the rest. The 
resulting function rh{p) is multivalued (since the theta functions are), but this is of 
no relevance for our purpose here. 

Proof. We will drop the dependence on t for notational simplicity. By p.9p the 
condition m{pi) = (O 0) implies that the Jost solutions 4'-{zi,n) and '0-f(zi,n) 
are linearly dependent. This can only happen, at a band edge, zi G {Ej}, or at an 
eigenvalue zi = pj. 

We begin with the case zi — pj. In this case the derivative of the Wron- 
skian W{z) — a{n)(ip^{z,n)ip-{z,n + 1) — ip^{z,n+ l)ip^{z,n)) does not vanish 
■j^W{z)\z=zi ^ (^ (6.11)]). Moreover, the diagonal Green's function g[z,n) — 
W{z)~^'ip+{z,n)ip-{z,n) is Herglotz and hence can have at most a simple zero at 
z = Pj. Hence, if ip^lpj^n) = ^^{pj,n) = 0, both can have at most a simple zero 
at z = Pj. But T{z) has a simple pole at pj and hence T{z)'ip^{z,n) cannot vanish 
at z — Pj, a contradiction. 

It remains to show that one zero is simple in the case zi G {Ej}. To show this 
we will introduce a local coordinate C = \/z — zi and show that ^M^(C)lc=o ^ in 
this case as follows: We will consider every function of 2; as a function of C and, by 
abuse of notation, denote it by the same name. Moreover, we will assume zi 7^ pj 
for all j. Then, note that V-'±(C) (where / denotes the derivative with respect to 
C) again solves H^'^{0) — zi-0j_(O) if zi G {Ej}. Moreover, by W{0) = we have 
'0+(O) — cijj-lO) for some constant c (independent of n). Thus we can compute 

w^'(o) = w^(^;(o),^_(o)) + i^(0+(o),^L(o)) 

= c-1i^(7A;(O),^+(O)) + cM^(0_(O),V/ (0)) 



14 



H. KRUGER, AND G. TESCHL 



by letting n — > +00 for the first and n — > —00 for the second Wronskian (in which 
case we can replace V'±(0,n) by V'g,±(0, n)), which gives (cf. [HI Chap. 6]) 

W (0) = -rpfo — hm — Tt 0. 

2nf=i(^i-M,)C-o C 

Hence the Wronskian has a simple zero. But if both functions had more than simple 
zeros, so would the Wronskian, a contradiction. The case where zi — fij for some j 
is similar. Just observe that in this case both 'ip± as well as ipq,± have singularities, 
which however cancel in the quotient appearing in rh{p). □ 

4. A UNIQUENESS RESULT FOR SYMMETRIC VECTOR RiEMANN-HlLBERT 

PROBLEMS 

In this section we want to investigate uniqueness for the meromorphic vector 
Riemann-Hilbert problem 

(7711) > -P^-, {m2)>~V^ 

1' 



(4.1) m{p*) — ■m{p) 



1 



m{co+) = (1 7712(00+)) . 

where S is a nice oriented contour (see Hvpothesis lA.ip . symmetric with respect to 
p 1-^ p* , and V is continuous satisfying 

(4.2) f(p*)=(j J)' pes. 

The normalization used here will be more convenient than (|3.23p . In fact, p.23p 
will be satisfied by 7773 {oo^)m{p). 

Now we are ready to show that the symmetry condition in fact guarantees unique- 
ness. 

Theorem 4.1. Suppose there exists a solution m{p) of the Riemann-Hilbert prob- 
lem (|4.ip such that for 

'9{z{p*j) 

e{z{p))^ 

the equality rh{p) = (O O) can happen at most for p e {^-j} in which case 

limsup^^^;^. j-^^ is bounded from any direction for j = 1 or j = 2. 

Then the Riemann-Hilbert problem (|4.ip with norming condition replaced by 

(4.4) m{oo+) = (a 7772(00+)) 

for given a G C, has a unique solution ma{z) — am{z). 

Again, the theta functions in (|4.3p are used to cancel the poles at fi and fi* such 
that their contribution is separated from the rest. 

Proof. Let ma(z) be a solution of (|4.ip normalized according to l|4.4p . Then we 
can construct a matrix valued solution via M = (777, 777c,) and there are two possible 
cases: Either det M(z) is nonzero for some z or it vanishes identically. 



(4.3) rh{p) = m{p) 
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We start with the first case. By the Lemma 13.41 we can rewrite ah poles as 
jumps with determinant one. Hence, the determinant of this modified Riemann- 
Hilbert problem has no jump. Hence it is a meromorphic function whose divisor 
satisfies (det(A/)) > — — I^/i*. Since T)jx is a nonspecial divisor, so is Vjx + Vj^* 
and the Riemann-Roch theorem implies that there are g + 1 linearly independent 
meromorphic functions of this kind. By inspection they are given by 

(4.5) det(M(p)) = , z = Trip), 

where P{z) is a polynomial of degree at most g. But taking determinants in 

gives a contradiction. 

It remains to investigate the case where det(Af) = 0. In this case we have 
uiaip) = 5{p)m{p) with a scalar function 5. Moreover, 5{p) must be holomorphic 
for z G C\I] and continuous for z G S except possibly at the points where m(zo) = 
(O O). Since it has no jump across E, 

6+{p)m+{p) ^ nia^+ip) = ma,^{p)v{p) = S-{p)m^{p)v{p) = 6-{p)m+{p), 

it is meromorphic with divisor 

3 

But the symmetry 5{p) — 5{p*) requires at least second order poles at the branch 
points (if any at all), which shows that 6 is constant. This finishes the proof. □ 

Furthermore, note that the requirements cannot be relaxed to allow (e.g.) second 
order zeros in stead of simple zeros. In fact, if m{p) is a solution for which both 
components vanish of second order at, say, p = Ej, then rfi(p) = ^J^. m{p) is a 
nontrivial symmetric solution of the vanishing problem (i.e. for a = 0). 

By Lemma 13.61 we have 

Corollary 4.2. The function rn{p, n, t) defined in (13. 9p is the only solution of the 
vector Riemann-Hilbert problem (j3.19p - (j3.23p . 



5. The stationary phase points and the nonlinear dispersion relation 

In this section we want to look at the relation between the energy A of the 
underlying Lax operator Hq and the propagation speed at which the corresponding 
parts of the Toda lattice travel, that is, the analog of the classical dispersion relation. 
If we set 



(5.1) v{X) ^ lim 



(A+ie, + ) 



the nonlinear dispersion relation is given by 
(5.2) viX) = J. 

Recall that the Abelian differentials are given by p.l4p and (|2.15p . 
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For p £ R\a{Hq) we have 
(5.3) v{p) 
that is, 



r(p,+)o 



U!, 



(5.4) = 7 ^ Hp,j)^0. 

In other words, v{p) is precisely the velocity of a soliton corresponding to the 
eigenvalue p. 

For A € (^{Hq) both nominator and denominator vanish on ij{Hq). Hence by de 
I'Hospital we get 

that is, 

(5.6) = 7 ^ 0'(A,y)=O, 

where 0, defined in (|3.14p . is the phase of factorization problem (|3.13p . In other 
words, v{X) = J if and only if A is a stationary phase point. 

Invoking (|2.14p and (|2.15p , we see that the stationary phase points are given by 

(5.7) fi(^-^^)+7n(^-^^)=o- 

Due to the normalization of our Abelian differentials, the numbers Aj, I < j < g, 
are real and different with precisely one lying in each spectral gap, say Xj in the 
j'th gap. Similarly, Aj, < j < are real and different and Xj, 1 < j < g, sits in 
the j'th gap. However Aq can be anywhere (see [211 Sect. 13.5]). 

Lemma 5.1 ([141). Denote by Zj{v), < j < g, the stationary phase points, where 
V — J. Set Aq = — oo and Xg+i ~ oo, then 

(5.8) Xj < zj{v) < Xj+i 

and there is always at least one stationary phase point in the j 'th spectral gap. 
Moreover, Zj{v) is monotone decreasing with 

(5.9) lim Zj{v) = A,+i and lim Zj(v) = Xj. 

V — * — oo V — ^oo 

So, depending onn/t there is at most one single stationary phase point belonging 
to the union of the bands a{Hq), say Zj{v). 
We now can establish that ^(A) is monotone. 

Lemma 5.2. The function v{p) defined in (|5.ip is continuous and strictly mono- 
tone decreasing. Moreover, it is a bijection from M to M. 

Proof. First of all observe that w(A) is continuous. This is obvious except at the 
band edges X = Ej. However, computing limA^B^ v{X) using again de I'Hospital 
establishes continuity at these points as well. 
Furthermore, for large p we have 

(5.10) lim 



pHoo -p/l0g(|p|) 
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which shows hmp^±oo v{p) — =poo. 

In the regions, where there is one stationary phase point Zj{v) G <^{Hq) we know 
that Zj{v) is the inverse of v{X) and monotonicity follows from the previous lemmas. 
In the other regions we obtain from v{C,{z)) = z by the implicit function theorem 
(5.11) 

^/__pl/2 ^^}Eo ^°o+°0- -_??l/2 ^^x ^^Jgp ^°g+°0- 



which shows strict monotonicity since Re /^'^'+Loo+oo_ > for C e ^ViHq) and 



^jlf) < C('^) < ^j-i(^) for C(^^) € {E2j-i, E2j) (if we set z_i = oo, Zg+i — —oo, 
E-i — — oo, E2g+2 = oo). To see the last claim we can argue as follows: If ({v) were 
below Zj{v) at some point it would decrease as v decreases whereas Zj{v) increases 
as V decreases. This contradicts the fact that both must hit at i?2j-i. Similarly we 
see that ({v) stays below Zj{v). □ 

In summary, we can define a function ((n/t) via 

(5.12) v{C) = ^. 

In particular, different solitons travel at different speeds and don't collide with each 
other or the parts corresponding to the continuous spectrum. 

Moreover, there is some Co for which w(Co) = and hence there can be stationary 
solitons provided Co ^ '^{Hq)- To show that this can indeed happen and to shed 
some further light on the location of this point we establish the following facts: 

Lemma 5.3. (i) satisfies 



(5.13) _ _|cr(i7g)| < Ao < — +21^(^9) 



where \a{Hq)\ = X]f=o(^2j+i ~ -^2j) denotes the Lebesgue measure of 
a{H,). 

(ii) There exists a unique Co such that v{Co) — which lies inside the convex 
hull of the spectrum a{Hq). Furthermore, if Co G cr(-ffo) or Aq G a{Ha), 
then Co = Aq. 

Proof, (i). Observe that by (|TT5| 

J=0 j = l j=l ^ 

and the claim follows using A^ e (-£'2j_i, £'2j)j j = ■ ■ ■ t 9- 

(ii). Existence and uniqueness of Co follows since u is a bijection. To check that Co 
is in the convex hull of a{Hq), it suffices to check that v{Eo) > and v{E2g+i) < 0, 
this follows from (ISJ]) using Xj,Xj G (£^2^-1, i?2j), and (I5.13p . That Co = Aq if 
C € cr{Hq) or Aq e cr{Ho) follows again from (|5.7p and Xj G (£'2^-1, i?2j). □ 

This lemma implies in particular that Co G ^{Hq) if and only if Ao G a{Hq). To 
see that both possibilities can occur, observe that if the band edges are symmetric 
with respect to A 1-^ —A, the same will be true for Xj and Xj. 
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6. The partial transmission coefficient 
Define a divisor V^^n t) of degree g via 

(6-1) aEg{'D^(n,t)) = aEo('^^(n,t))+S{n/t), 

wfiere 

(6.2) Se{n/t) = 2 ^ A^^i^p,) + ± f \ogil ~ \R\')Q, 

P.<an/t) 27ri7e(„/*) 

wfiere C{n/t) = S n 7r~i((-oo, C(n/0) and Cin/t) is defined in ([51^ . 

Tfien 'Di^(^n,t) is nonspecial and T:{Oj{n,t)) — Vj(n,t) € M with preciseiy one in 
each spectrai gap (see [l4]). 

Then we define the partiai transmission coefficient as 
(6.3) 

,'e(^(n,t)+^(nA))0(z(n-l,<)+^(nA))V/' 9iz{p,n,t)) 
T{p,n,t) = > 



0{z{n,t)) e{z{n~l,t)) J 0{zip,n,t)+S{n/t)) 

TT exp ( - [ ujp^pi) I exp ( [ iog(l - \R\^)ujpEo ] , 

where ^(n, t) is defined in (|6.2p and ujpq is the Abeiian differential of the third Icind 
with poles at p and q. In the case where we have the full transmission coefficient 
this formula was derived in [23] (see also [5] and [7]). 

The function T{p,n,t) is meromorphic in M \ C(n/t) with first order poles at 
Pk < Q{n/t), i)j{n,t) and first order zeros at fij{n,t). 

Lemma 6.1. T{p,n,t) satisfies the following scalar meromorphic Riemann-Hilbert 
problem: 

T+{p,n,t)=T^{p,n,t)il~\Rip)f), peC{n/t), 
(6.4) iT{p,n,t))= Vp.~ X>p, +X>^(„_t) -X>^(„_4), 

Pfc<C(n/t) Pk<C{n/t) 

r(oo-|_, n, t)r(cx3_, rt, i) — 1, T{GO+,n,t) > 0. 

Moreover, 
(i) 

g _ 

Tip*,n,t)Tip,n,t)^Y[^-^, z = nip). 



(ii) T{p, n, t) ~ T{p, n, t) and in particular T(p, n, t) is real-valued for p € T,. 
Proof. This can be shown as in [14, Thm. 4.3]. □ 

We will also need the expansion around oo+ given by 

(6.5) T{p,n,t)=Toin,t)fl + ^^^^ + Oi^)\, P={z,+), 
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where 



To(n,<) =r(oo+,n, i) = 



6>(z(n, t)) 9{z{n - l,t) + 5{n/t)) \ ^'"^ 
e{z{n,t)+5{n/t)) 6{z{n~l,t)) 



(6.6) 



and 



1 d ^ s) +5{n,t)) 



2ds 



e{z{n,s)) 



where VLq is the Abehan differential of the second kind defined in (|2.f4p . This 
follows as in [23l Sect. 4]. 



7. SOLITONS AND THE SOLITON REGION 

This section demonstrates the basic method of passing from a Riemann-Hilbert 
problem involving solitons to one without. Solitons are represented in a Riemann- 
Hilbert problem by pole conditions, for this reason we will further study how poles 
can be dealt with in this section. We follow closely the presentation in Section 4 of 

In order to remove the poles there are two cases to distinguish. If pj > C,{n/t) 
the jump is exponentially close to the identity and there is nothing to do. 

Otherwise we need to use conjugation to turn the jumps into this form exponen- 
tially decaying ones, again following Deift, Kamvissis, Kriecherbauer, and Zhou [4]. 
It turns out that we will have to handle the poles at pj and p* in one step in order 
to preserve symmetry and in order to not add additional poles elsewhere. 

For easy reference we note the following result which can be checked by a straight- 
forward calculation. 

Lemma 7.1 (Conjugation). Assume that ECS. Let D he a matrix of the form 

<") -<^) - {V X,) • 

where d : M\S C is a sectionally analytic function. Set 

(7.2) m(p) ~ m{p)D{p), 
then the jump matrix transforms according to 

(7.3) v{p) ^ D_{p)-^v{p)D+{p). 

rh{p) satisfies p.22p if and only if m{p) does. Furthermore, rh(p) satisfies p.23p . 
if m{p) satisfies p.23|) and d{oo±)d(oozf) = 1. 

In contradistinction to [TB], we will no longer have det{D{p)) — 1, but det(-D(p)) 
will be the lift of a rational function with g zeros and poles. 
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Lemma 7.2. Introduce 

(7.4) B{p,p)^C,{n,t) , .. Bip,p). 

e{z{p,n,t) + 2Aj^^{p)) 

Then B(.,p) is a well defined meromorphic function, with divisor 

(7.5) iB{.,p)) = -V,+V^-Vp. +Vp, 
where v is defined via 

(7.6) {V,) = (2?^ ) + '^Ae, (p) ■ 
Furthermore, B{.,p) has a pole at p* , and 

(7.7) S(oo+,p)B((X3_,p) = l, 
^f 

(7.8) Cp(n,t) ^ eiz{n,t)) 0(z(n-l,<)) ' 

Proof. Wc start by checking single valuedness. The a-periods follow from normal- 
ization. For the b periods, we compute for E E {E21-1, E21) 

B{E + ie,p) 
TiQ B{E - ie, p) 

using (A. 21) in 21j. Now using (A. 68) in [21 the claim follows. 

The normalization condition (|7.7p follows by a computation using (|2.17p . □ 

Now, we can show how to conjugate the jump corresponding to one eigenvalue. 

Lemma 7.3. Assume that the Riemann-Hilbert problem for m has jump conditions 
near p and p* given by 



% ' = exp (27ri(2ABo,;(p))) 



(7.9) 



/ 1 0\ 

m+{p) = TO_(p) 7(p) ^ \ , pe ^e{p), 

( Y lip*) \ 



m+(p)=TO_(p)l Av)-P\ peE,(p*), 



and satisfies a divisor condition 

(7.10) (mi)>-I?^., {m2)>-Vf^. 

Then this Riemann-Hilbert problem is equivalent to a Riemann-Hilbert problem for 
rh which has jump conditions near p and p* given by 

A B(p,p-)(7r(p)-p) \ 

m+(p) — rh_{p) j{p)B(p-',p'-) 1 p£ T,^{p), 

(7.11) ; ^ ' . 

/ 1 o\ 

m+{p) ^ m-(p) Bip',p')i^(p)-p) , , peT,,{p*), 

\ l{p')B(p,p') J 

divisor condition 

(7.12) {mi)>-Vo*, {m2)>-Vo, 
where I?j> is defined via 

(7.13) aE^,(Po) - a^„(I?p) + 2Ae^{p). 
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and all remaining data conjugated (as in Lemma pOp by 



Proof. Denote by U the interior of (p) . To turn 7 into 7 




Dip) 



and note that D{p) is meromorphic away from the two circles. Now set rh{p) = 
m(jp)D{p). The claim about the divisors follows from noting, where the poles of 
B{p,p) are. □ 

This result can be applied iteratively to conjugate all eigenvalues pj < Q{n/t) 
as follows: One starts with the original poles p = pP and applies the lemma with 
p — Pi resulting in new poles p^ = v. Then one repeats this with p = p^, p = P2, 
and so on. Finally, we will also need a last conjugation step to factor our jump 
matrices into upper and lower triangular parts as demonstrated in [T3] . Combining 
all steps we end up with the following conjugation: 

Abbreviate 



and introduce 



D{p) 



lk{p,n,t) 



1 

7fc(P:",*) 

7r(p)-Pfc 





R 



1/2 
23+2 



(Pk) ipq{p,n,t) 
ULiiPk - Pi) '4'q{P*,n,t) 

'{p)-pk 



Ik 



\'rk{p',n,t) 




where 



Doip) 



T{p\n,t) 



pen+ 

k(p)-Pfcl<e 

pen_ ' 
else, 







Pk < Cin/t), 



r(p,n,i), 
Note that D{p) is meromorphic in M\I](zo) and that we have 



D{p* 



1 

1 



Dip) 



1 

1 



Now we conjugate our problem using D(p) and observe that, since T{p,n,t) has 
the same behaviour as T{p) for p a band edge, the new vector rh{p) = m{p)D{p) is 
again continuous near the band edges. 
Then, the divisor conditions are shifted 



(7.15) 



("T-i) > -2?s(„,t)*' ("^-2) > 
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Moreover, using Lemma [7. II and Lemma [7.31 the jmiip corresponding pk < C,{n/t) 
(if any) is given by 

A T{p,n,t)(-K{p)-pk) \ 

(7.16) ) ^ 

'"{P) = T{p\n,t){^{p)-pk) . , pe Se(p^), 

\ -tk{p' ,n,t)T{p,n,t) ^) 

and corresponding > C("-/^) (if ^^^y) by 

, > \T(p,n,t)(7r(p)-pfe) 

/i 7fc(p%",*)T(p,n,t) \ 

V{p) = I J T(p*,n,tKT(p)-Pfe) J , p g 5]e(pD- 

In particular, all jumps corresponding to poles, except for possibly one if pk — 
({n/t), are exponentially decreasing. In this case we will keep the pole condition 
which now reads 

f~ , , , jk{p,n,t)T{p*,n,t) . , ,\ ^ „ 

[mi(p) + — 7m2[p) > -Vc(n.tr, near p^, 

(^y^g^ ^ T{p,n,t){TT{p) - Pk) I ' ' 

( lk{p*,n,t)T{p,n,t) ~ / ^ , ~ , A ^ ^ 

n, t) in{p) -p^r ^^P^^ ^ '^^("-*) ' ^^"^^^^ 

Furthermore, the jump along S is given by 



(7.19) v{p) = 

where 



L(p)-i6+(p), ^(p)>C(nA), 
B_{p)-^B+{p), n{p)<an/t), 



I , s A ^i?(p*)e(p*)e-*'^(p)\ r 



,0 1 i'"+^^^-i^i?(p)e(rte*^(^) 1 



and 



^-^^^^(-^i?(.)e(p)e*^(^) ij'^-^^^^-l^o 1 j 

In our next step we make a contour deformation and move the corresponding parts 
into regions where the off-diagonal terms are exponentially decreasing. For this we 
take some sufficiently small loops Ck around each spectral band [E2k,E2k+i] on the 
upper sheet (see Figure [1]). In particular, these loops must not intersect with any 
of the loops around the eigenvalues pj. Then, an investigation of the sign of Re(0) 
shows that 



(7.20) 



fRe((/.(p)) < 0, peDk, n{p) > ((n/t) 
[Re((/.(p)) > 0, pe Dk, nip) < dn/t) 
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Figure 1 . The small lens contour around a spectral band. Views 
from the top and bottom sheet. 



and we can deform our contour according to 



(7.21) 



m{p) 



= < 



rh{p)h+(pY 
fh{p)b-{p)- 
■m{p)B+{p) 
fh{v)B- (p) 
m{p), 



peDk, 7t{p) > C{n/t), 

peDl tt{p) > C{n/t), 

peDk, Trip) < C{n/t), 

peDl, tt{p) <C{n/t), 
else, 



such that the jump on S disappears and on lJfe=o(^'= ^ ^k) given by 



(7.22) 



v{p) = < 



b+ip), 

B+{p), 

vip), 



peCk, n{p) > Cin/t), 
peCl7T{p)>C{n/t), 
p£Ck, n{p) < Cin/t), 
p£Cl, TT{p)<C{n/t), 
else. 



The jumps on the small circles around the eigenvalues remain unchanged. 

Here we have assumed that R has an analytic extension to the corresponding 
regions. We will show how this restriction can be overcome below. Moreover, in 
order to obtain uniform errors we need to require that the Vj retain some minimal 
distance from the jump contour. But this implies that some regions of (n, t) values 
are not covered. However, we can evade this obstacle by slightly deforming our 
contour such that it has some positive distance to the first one. This way we also 
cover the missing regions. 

Now we are ready to prove Theorem II. II bv applying Theorem lA. 61 as follows: 

If \(^{n/t) — pk\ > £ for all k we can choose 70 = and Wq by removing all jumps 
corresponding to poles from w*. In particular, the error between the solutions of 

and Wq is exponentially small. This proves the second part of Theorem II . II upon 
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comparing 
(7.23) 

with ((37T51) using (|6.5p . 

Otherwise, if ICC*^/*) — Pfcl < ^ for some fc, we choose 7o = 7fc(n,t) and Wq = 0. 
Again we conclude that the error between the solutions of w* and Wq is exponentially 
small. 

If R{p) has no analytic extension, we will approximate R{p) by analytic functions 
in the spirit of [3]. In fact, as in [ITl Sect. 5] one sees that it indeed suffices to find 
an analytic approximation for the left and right reflection coefficients. Moreover, 
for each spectral band (viewed as a circle on the Riemann surface) one can take 
the imaginary part of the phase as a coordinate transform and then use the usual 
Fourier transform with respect to this coordinate (compare [171 Lem. 5.3]). In 
order to avoid problems when one of the poles Vj hits E, one just has to make the 
approximation in such a way that the nonanalytic rest vanishes at the band edges 
(cf. Remark I A. 7p . That is, split R according to 

Rip)=RiE,,) '~\ +RiE,,^,) 



E2j + l — E2j E2j — E2j + l 



(7.24) ±^z-E2,^z-E2j+iR{v), P=(^,±), 

and approximate R. 

Remark 7.4. Note that one can even do slightly better by using the weighted mea- 
sure — ii?2g+2(p)'^^ ''^ *^ which case it suffices if R is just C^^^'^^CE) for some 
7 > rather than C'+^(I]). In fact, one can show that the Cauchy operators are 
still bounded in this weighted Hilbert space (cf. pTl Thm. 4.1]j. 

Acknowledgments. G.T. gratefully acknowledges the extraordinary hospitality 
of the Department of Mathematics at Rice University, where part of this research 
was done. 

Appendix A. Singular integral equations 

In this section we show how to transform a meromorphic vector Riemann-Hilbert 
problem with simple poles at p, p*, 

m-^-{p) = m^{p)v{p), p G S, 

(mi) > -Vf,^ ~ Vp, (ma) > - Vp, , 

■^2g+2(p) 7j V'g(p) 



(A.l) (mi(p) - T^":. „r;i:> 2(p)) > -Vl near p. 



-^29^2 (P) 7 V'g(p) 

nLl(P - Mfc) 7r(p) - ptpq{p*)' 



mi{p) +m2{p)) > -Vp, near p* , 



m{p*) = m{p) ^ 
m(cx)+) = (l ma) , 



into a singular integral equation. Since we require the symmetry condition p.22p 
for our Riemann-Hilbert problems, we need to adapt the usual Cauchy kernel to 
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preserve this symmetry. Moreover, we keep the single sohton as an inhomogeneous 
term which wiU play the role of the leading asymptotics in our applications. 

Hypothesis H. A.l. Let S consist of a finite number of smooth oriented finite 
curves in M which intersect at most finitely many times with all intersections being 
transversal. Assume that the contour S does not contain oo± and is invariant 
under p ^ p* . It is oriented such that under the mapping p i—^ p* sequences 
converging from the positive sided to E are mapped to sequences converging to the 
negative side. Moreover, suppose the jump matrix v can be factorized according to 
V = bZ^bjf. — {I — w^)^^{I + w^), where w± — ±(6± — I) are continuous and satisfy 

(A.2) «^±b*)=(j ^eE. 

In order to respect the symmetry condition we will restrict our attention to the 
set Ls(E) of square integrable functions / : E — > such that 

(A.3) /(,)=/(.•)(; ; 

Clearly this will only be possible if we require our jump data to be symmetric as 
well (i.e., Hypothesis I A. II holds) . 

We begin by introducing the Cauchy operator following Section 5 in [T4| . Given 
a nonspecial divisor I?^, we introduce the Cauchy kernel 

(A.4) n§-''^c.,, + j24''ipKj, 

where 



(A.5) /f''(p) = ^c,,(A) f'i 

1=1 •'P 

Here q is the (normalized) Abelian differential of the second kind with a second 

order pole at q. Note that Ij''{p) has first order poles at the points fi. 
Introduce 

M.P _ / 



(A.6) 

and define the Cauchy operator by 



\ 4^% 



1 



(A.7) {Cf){p) = — / fisnl;'" 

acting on vector- valued functions / : E — > C^. We will assume that I?^ does not 
hit E (see Remark I A. 71 below for the case where this assumptions does not hold). 

Lemma A.2. Assume Hypothesis \A.l[ The Cauchy operator C has the properties, 
that the non-tangential boundary limits 

(A.8) {C±f){q)= hm / fnf 
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from the left and right of S respectively ( with respect to its orientation ) are hounded 
operators L'j.iY,) — > The bound can he chosen independent of the divisor as 

long as it stays some finite distance away from E. The operators C± satisfy 

(A.9) C+-C-^I 

and 

(A.IO) (C/)(p*)-(0 *), = 0). 

Furthermore, C restricts to Lg(I]), that is 

(A.ll) (C/)(p*) = (C/)(p)('? ]\, peM\S 



for f G and if w± satisfy (|A.2[) we also have 

(A.12) C±ifw^){p*) ^ C^ifw±)ip) (^^ J), pes. 

Proof. Follows from the properties of Cauchy operators, see Theorem 5.1 in [14] . □ 

We have thus obtained a Cauchy transform with the required properties. Fol- 
lowing Sections 7 and 8 of [1], we can solve our Riemann-Hilbert problem using 
this Cauchy operator. 

Introduce the operator : LKT,) — > Ll{Y^) by 

(A.13) C^f^C+{fw^) + C-{fw+), /eL2(S). 

Recall from Lemma 13.51 that the unique solution corresponding to w = I is given by 

mo(p) = (/(p*) /(p)), 
for some given /(p) with (/) > — P/i — Dp* . Since we assumed T>^ to be away from 
E, we clearly have niQ e L^(S). 

Theorem A. 3. Assume Hvvothesis \A.l[ 

Suppose m solves the Riemann-Hilbert problem (|A.1[) . Then 

(A.14) m(p) = {1 - co)mo{p) + J^fi{s){w+{s) + w-{s))i$''' , 

where 

/i = 771+67^ = m_6l^ and cq — ( / fj,{s){w+{s) + W-{s))fl^'\_] . 

Here {m)j denotes the j 'th component of a vector. Furthermore, /i solves 
(A.15) (I-C^)Ai= (l-co)mo(p). 

Conversely, suppose p. solves 
(A. 16) (I-C„)/i = mo, 

and 

Co = ^ fl{s){w+{s) + w-(s))fl^'j_^ 7^ 0, 

then m defined via (|A.14|) . (1 — cq) = (1 — cq)^^ and /i = (1 — cq)^"'^//, solves 
the Riemann-Hilbert problem (jA.ip aric? /i = m±b^ . 

Proof. Follows as in [TB]. □ 
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Hence we have a formula for the solution of our Ricmann-Hilbert problem m{z) 
in terms of (I— Cu,)~^too and this clearly raises the question of bounded invertibility 
of I — Cw This follows from Fredholm theory (cf. e.g. [25]): 

Lemma A. 4. Assume Hypothesis \ A. 1[ 

The operator I — Cw is Fredholm of index zero, 

(A.17) ind(I-C^) = 0. 

Proof. Follows as in [H], [16]. □ 

By the Fredholm alternative, it follows that to show the bounded invertibility of 
I — Cw we only need to show that ker(I — C^) = 0. The latter being equivalent to 
unique solvability of the corresponding vanishing Ricmann-Hilbert problem. 

Corollary A. 5. Assume Hypothesis \A.l\ 

A unique solution of the Riemann-Hilhert problem (jA.ip exists if and only if the 
corresponding vanishing Riemann-Hilhert problem, where the normalization condi- 
tion is replaced by to(0) — (O rn2) , has at most one solution. 

We are interested in comparing two Riemann-Hilbert problems associated with 
respective jumps wq and w with \\w — wo||oo small, where 

(A. 18) ll^lloo = ||w+||l~(s) + ||w_||loo(s). 

For such a situation we have the following result: 
Theorem A. 6. Assume that for some data Wq the operator 
(A.19) I-C^. :L2(£)^L2(E) 

has a bounded inverse, where the bound is independent oft, and let C = Co; 7* = 7o- 
Furthermore, assume w* satisfies 

(A.20) \\w* - wlW^ < ait) 

for some function a{t) —> as t oo. Then (I — C^t)"^ : Ls(S) ^s(^) '^^^'^ 
exists for sufficiently large t and the associated solutions of the Riemann-Hilbert 
problems (jA.ip only differ by 0{a{i)). 

Proof. Follows as in [TB]. □ 

Remark A. 7. The case where one (or more) of the poles fij lies on S can be 
included if one assumes that w± has a first order zero at fij. In fact, in this case 
one can replace ii{s) by /i(s) = (7r(s) — iJ,j)ii{s) and w±(s) by 'w±{s) — (tt{s) — 
fj.j)-^w±{s). 

Acknowledgments. We thank A. Mikikits-Leitner for pointing out errors in a 
previous version of this article. G.T. gratefully acknowledges the extraordinary 
hospitality of the Department of Mathematics at Rice University, where part of 
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